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We give detailed description of the transport spin current in the chiral helimagnet. Under the 
static magnetic field applied perpendicular to the helical axis, the magnetic kink crystal (chiral 
soliton lattice) is formed. Once the kink crystal begins to move under the Galilean boost, the 
spin-density accumulation occurs inside each kink and there emerges periodic arrays of the induced 
magnetic dipoles carrying the transport spin current. The coherent motion of the kink crystal 
dynamically generates the spontaneous demagnetization field. This mechanism is analogous to the 
Doring-Becker-Kittel mechanism of the domain wall motion in ferromagnets. To describe the kink 
crystal motion, we took account of not only the tangential (^-fluctuations but the longitudinal 6- 
fluctuations around the helimagnetic configuration. Based on the collective coordinate method and 
the Dirac's canonical formulation for the singular Lagrangian system, we derived the closed formulae 
for the mass, spin current and induced magnetic dipole moment accompanied with the kink crystal 
motion. To materialize the theoretical model presented here, symmetry-adapted material synthesis 
would be required, where the interplay of crystallographic and magnetic chirality plays a key role 
there. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The core problem in the multidisciplinary field of 
spintronics is how to create, transport, and manipu- 
late spin currents^ The key notions include the current- 
driven spin-transfer torque 2 ^^^ and resultant force 
acting on a domain wall (DW)£i£ in metallic fer- 
romagnetic/nonmagnetic multilayers, the dissipation- 
less spin currents in paramagnetic spin-orbit coupled 
systems ^ 10 ' 11 and magnon transport in textured mag- 
netic structures^ A fundamental query behind the issue 
is how to describe transport spin currents^ To make 
clear the meaning of the spin currents, we need to note 
the spin can appear in the macroscopic Maxwell equa- 
tions only in the form of spin magnetization. In this 
viewpoint, the spin current is understood as the devi- 
ation of the spin projection from its equilibrium value. 
An emergence of the coherent collective transport in non- 
equilibrium state is then a manifestation of the dynamical 
off-diagonal long range order (DODLRO) j 14 ' 15 

On the other hand, the physical currents are classified 
into two categories, i.e., the gauge current originating 
from the gauge invariance and the incrtial current origi- 
nating from the Galilean invariance. The electric current 
is the gauge current, where the electric charge is coupled 
to the electromagnetic U(l) gauge field. The electro- 
magnetic field is a physical gauge field that has its own 
dynamics, i.e., we know the electromagnetic field energy. 
Then, the charge current ji and the charge density p are 
related via the continuity equation dp/dt = —dji/dxi. 
On the other hand, a typical example of the inertial cur- 
rent is the momentum current in a classical ideal fluid, 
where the momentum current Ily satisfies the continu- 
ity equation, d(pvi)/dt = —dUij/dxj, and given by 
Tlij — PSij + pviVj with P being equilibrium pressured 



The non-equilibrium current is described by pviVj. In 
the spin current problem, at present, we have no known 
gauge field directly coupled to the spin current. There- 
fore, a promising candidate is the inertial current of the 
magnetization. 

Historically, Doringil pointed out that the longitudi- 
nal component of the slanted magnetic moment inside 
the Bloch DW emerges as a consequence of translational 
motion of the DW. An additional magnetic energy asso- 
ciated with the resultant demagnetization field is inter- 
preted as the kinetic energy of the wall. This idea was 
simplified by Becker— and Kittel^ Recent progress of 
material synthesis sheds new light on this problem. In a 
series of magnets belonging to chiral space group with- 
out any rotoinversion symmetry elements, the crystallo- 
graphic chirality gives rise to the asymmetric Dzyaloshin- 
skii interaction that stabilizes either left-handed or right- 
handed chiral magnetic structures^ In these chiral he- 
limagnets, magnetic field applied perpendicular to the 
helical axis stabilizes a periodic array of DWs with def- 
inite spin chirality forming kink crystal or chiral soliton 
latticed 

We recently proposed a new way to generate a spin cur- 
rent in the chiral hclimagncts with magnetic field applied 
in the plain of rotation of magnetization^ 2 - The mech- 
anism is quite analogous to the Doring-Becker-Kittel 
mechanism. We showed that the periodic spin accumu- 
lation occurs as a dynamical effect caused by the moving 
magnetic kink crystal (chiral soliton lattice) formed in 
the chiral helimagnet under the static magnetic field ap- 
plied perpendicular to the helical axis. The current is 
inertial flow triggered by the Galilean boost of the kink 
crystal. An emergence of the transport magnetic cur- 
rents is then a consequence of the dynamical off-diagonal 
long range order along the helical axis. 

In this paper, we give an extension of the results 
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touched on in Ref. 22 . In Sec. II, we give an overview 
of basic properties of the chiral magnets that materialize 
the theoretical model considered in this paper. In Sec. 
Ill, we present standard description of the kink crystal 
formation, and the vibrational modes around the kink 
crystal state. In Sec. IV, we apply the collective coordi- 
nate method to the moving kink crystal that makes clear 
the physical meaning of the mass and the magnon cur- 
rent carried by the moving system. In Sec. V, we perform 
quantitative estimates of the mass, magnetic current, and 
net magnetization induced by the movement. In Sec. VI, 
we discuss issues closely related to the present problem, 
including the background spin current problem, spin su- 
percurrent in the superfluid 3 He, and experimental as- 
pects of our effects. Finally, we summarize the paper in 
Sec. VII. 



II. CHIRAL HELIMAGNET 

In this section, we briefly review basic properties 
of chiral helimagnets that materialize our theoretical 
model. Recent progress of material synthesis promotes 
systematic researches on a series of magnets belonging 
to chiral space group without any rotoinversion sym- 
metry elements^ In the chiral magnets, the crystallo- 
graphic chirality possibly gives rise to the asymmetric 
Dzyaloshinskii interaction that stabilizes the chiral he- 
limagnetic structure, where either left-handed or right- 
handed magnetic chiral helix is formed. 20 As we will see, 
in the chiral helimagnets, magnetic held applied perpen- 
dicular to the helical axis stabilizes a periodic array of 
DWs with definite spin chirality forming kink crystal or 
chiral soliton latticed 

The chiral helimagnetic structure is an incommensu- 
rate magnetic structure with a single propagation vec- 
tor k = (0,0, A;). The space group Q consists of the 
elements {gi}. Among them, some elements leave the 
propagation vector ko = (0,0, k) invariant, i.e., these 
elements form the little group Gk i 23 i 24 The magnetic 
representation 2 ^ T mag is written as T rnag = r pcrm <g)r ax i al , 
where r pcrm and r ax i a i represent the Wyckoff permuta- 
tion representation and the axial vector representation, 
respectively. Then, T mag is decomposed into the non- 
zero irreducible representations of Gk • The incommen- 
surate magnetic structure is determined by a "magnetic 
basis frame" of an axial vector space and the propaga- 
tion vector k . In specific magnetic ion, the decom- 
position becomes T mag = ^^7?,^, where I\ is the ir- 
reducible representations of Gk . Then, we have two 
cases leading to the chiral helimagnetic magnetic struc- 
ture. Case I: The magnetic moments are described by 
two independent one-dimensional representations that 
form two-dimensional basis frames, or Case II: The mag- 
netic moments are described by a single two-dimensional 
representations that form two-dimensional basis frames. 
In these cases, the symmetry condition allows the chi- 
ral helimagnetic structure to be realized. Then, the 



structure is stabilized by the generalized Dzyaloshinskii 
interaction. The generalized Dzyaloshinskii interaction 
means symmetry-adapted anti-symmetric exchange in- 
teraction, not restricted to conventional Dzyaloshinskii- 
Moriya (DM) interaction caused by the on-site spin-orbit 
coupling and the inter-site exchange interactions. The 
presence of this term is justified by the existence of the 
Lifshitz invariant 2 ^ for the little group Gk ■ 

Among the inorganic chiral helimagnets, the best 
known example is the metallic helimagnet MnSi (T c ~ 
30K) that belongs to the cubic space group P2 1 3(a = 
4.558A)3& The metallic helimagnet Cr 1/3 NbS 2 (T c ~ 
120K) belongs to the hexagonal space group P6322 
(a = 5.75A, c = 12.12A)i22 The insulating copper 
metaborate, CUB2O4 (T c ~ 10K) has a larger unit 
cell and belongs to the tetragonal space group I42d 
(a = 11.48A, c = 5.620A)i2&2ii As examples of molecular- 
based magnets, the structurally characterized green nee- 
dle, [Cr(CN) 6 ][Mn(5 or i?)-pnH(H 2 0)]H 2 (T c ~ 38K), 
belongs to the orthorhombic space group P2\2\l\ (a = 
7.628A, b = 14.51A, c = 14.93A). The yellow needle, 
K .4[Cr(CN) 6 ][Mn(5)-pn](5)-pnH . 6 : {(S)-pn = (S)-l,2- 
diaminopropane) (T c ~ 53K), belongs to the hexagonal 
space group P61 (a = 14.77A, c = 17.57A}, 2 ! From the 
symmetry-based viewpoints, these space groups are all 
eligible to realize the chiral helimagnetic order. 



III. KINK CRYSTAL AND VIBRATIONAL 
MODES AROUND THE KINK-CRYSTAL STATE 

As shown in Fig. [TJ we consider a system of the chiral 
helimagnetic chains described by the model Hamiltonian 

n = -j s« ■ s s 

<i,3> 

+ D SiXSj-H-^Si, (1) 

<i,j> i 

where the first term represents the ferromagnetic cou- 
pling with the strength J > between the nearest neigh- 
bor spins Si and Sj. The second term represents the 
parity-violating Dzyaloshinskii interaction between the 
nearest neighbors, characterized by the the mono-axial 
vector D = De x along a certain crystallographic chiral 
axis (taken as the a;-axis). The third term is the Zee- 
man coupling with the magnetic field H = g/is_ffe y ap- 
plied perpendicular to the chiral axis. When we treat the 
model Hamiltonian ([1}, we implicitly assume that the 
magnetic atoms form a cubic lattice and the uniform fer- 
romagnetic coupling exists between the adjacent chains 
to stabilize the long-range order. Then, the Hamiltonian 
(TT]) is interpreted as a quasi one-dimensional model based 
on the interchain mean field pictured 

When H = 0, the long-period incommensurate heli- 
magnetic structure is stabilized with the definite chiral- 
ity (left-handed or right-handed) fixed by the direction 
of the mono-axial D-vector. The Hamiltonian (fTJ) is the 
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same as the model treated by Liu 3 -^ except that we ignore 
the single-ion anisotropy energy. Once we take into ac- 
count the easy-axis type anisotropy term, — K J2i(Sf ) 2 j 
the mean field ground state configuration becomes cither 
the chiral helimagnet for K < D 2 /J, or the Ising ferro- 
magnet for K > D 2 /J. In this paper, we assume K = 
and left an effect of K for a future study. 
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The magnetic kink crystal phase is described by the 
stationary soliton solution minimizing TL 1 9 = 7r/2 and 
cos [<fo(x)/2] = sn [mx/n, k], where sn is the Jacobi ellip- 
tic function with the elliptic modulus k (0 < k 2 < I). 25 i 32 
This solution corresponds to a periodic regular array of 
the magnetic kinks with the "topological charge" den- 

sity d x ipo (x) — 2 — dn I — x, k) as shown in Figs. ETb) 
k V k J 

and (c). The elliptic modulus k is found from the min- 
imization of energy per unit length that yields k/to = 
4E(K)/irqo£2. The period of the soliton lattice is given by 



^kink 



2kK(k) _ 8K(k)E(k) 



(3) 



FIG. 1: Schematic view of the model chiral helimagnet con- 
sidered here. 



Taking the semiclassical parametrization of 
Heisenberg spins in the continuum limit S (x) — 
S(cos9 (x) , sin 9 (x) cos(p(x) , sin 9 (x) simp (x)) by using 
the slowly varying polar angles 9{x) and ip(x) [see 
Fig. HJa)] , the Hamiltonian acquires the form 



H[<p{x),6(x)] 

rL 



= JS 2 



dx 



^{d x 9(x)} 2 + ^sin 2 9{d x <p(x)} 2 



-q sin 2 9 (x)d x ip(x) ~ m 2 sin9 (x)cos(p(x)] , (2) Jh/H c = k/E(k). 



where K (k) and E(k) denote the elliptic integrals of the 
first and second kind, respectively. The period increases 
from 2Tr/qo to infinity as k increases from zero to unity. 
In the limit of k, — > 0, the sn function approaches sin and 
/t/m— * 2/(70, i.e. ipo(x) — q x as it should be in the case 
of zero field. 

In the Hamiltonian @, the exchange processes favor 
the incommensurate (IC) chiral helimagnetic order, while 
the Zeeman term favors the commensurate (C) phase. 
The C-IC transition occurs at k = 1 provided E(l) = 
1, and the critical value of m is given by 7rq , o/4m c — 
1 i 33 i 34 i 35 The critical field strength H c is determined from 



where m = y H / JS, and L denotes the linear dimension 
of the system. From now on, all distances are measured 
in the lattice unit ao- The helical pitch in the zero field 
(to = 0) is given by qo = D/J. 
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Next, we consider the fluctuations around the kink 
crystal state. The studies of collective excitations in 
the system have been focused on the phasons (0-modc) 
presenting bending waves of domain walls of the soliton 
latticed In our analysis we are interested in the 0-modes 
also. We derive the spectrum of elementary exciations 
holding the scheme outlined in Refi^i 

In what follows, it is convenient to work with the di- 
mensionless coordinate 



2K («) 



4ink 4£ («) 



q a x. 



(4) 



We introduce L = mL/n and go 
the Hamiltonian ([2]) as 



Kqo/m, and rewrite 



FIG. 2: (a) Polar coordinates in the laboratory frame, (b) 
Formation of the magnetic kink crystal in the chiral helimag- 
nets under the transverse magnetic field, and (c) concomitant 
phase modulation. In (b), we depict a linear array of the spins 
along one chiral axis that is ferromagnetically coupled to the 
neighboring arrays. 
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2K (k) 

-^kink 



H, 



(5) 



where the dimensionless Hamiltonian TL is defined by 
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n = 



dx 



\{d s 6(x)f + ^sm 2 9(x) {d x ^{x)f 



qo sin 2 9 (x) d- <p (x) — k 2 sin 9 (x) cos <p (x) 



As the magnetic field increases from H = to H = 
H c , the parameter q~o = AE(k)/t: monotonously decreases 
from q~o =2 to q~o = 4/7r~1.273. The fluctuations consist 
of the vibrational (phonon) modes and the translational 
mode, that are separately treated. In this section, we 
examine the phonon modes. We write 

_ _ _ _ 

<p(x) = <po(x) +v(x), 6(x) = -+u(x) (7) 

and expand ([5]) up to u 2 and v 2 . Then we have Ji 

= J L dx(H +H u +H v +Hint)+O (u 2 , v 2 ) , where H cor- 
responds to the stationary solution. The interaction part 
contains —u 2 (d s v) 2 /2 and u A terms that are neglected 

here. The vibrational term V = f Q L dx (Ti u + H v ) is 
given by H u = uL u u, and H v — vL v v, where the dif- 
ferential operators, C v and C u , are defined by 

{Ai = + -K 2 COS^ > 

1 I I 

&• = - + cos<A) + 2 A ^ 2 ' 

The "gap function" reads as 

A(5) = ^2q {d s tpo) - (d s ip ) 2 

= 2\j qodn (x, k) — dn 2 (x, n), (9) 

where the relation dx<po = 2dn (x, k) was used. The 
minimum and maximum value of the gap are given by 

A max = go, A min = A(K) = 2V^^q - «/, (10) 

respectively, where k' = yl — k 2 is the complemen- 
tary modulus. We see that the gap closes at the C-IC 
transition. In Figj3] (a), we show the spatial variation 
of the gap function. The n dependence of the mini- 
mum gap is shown in Fig[3] (b). For small k, we have 
A max ~ 2-k 2 /2-3k 4 /32, and A min ~ 2-k 2 /2-7k 4 /32. 
Therefore, A max /A min ~ 1 and it is appropriate to ap- 
proximate A(x) ~ 2 for the case of weak field. This 
approximation amounts to approximating dn (x, k) ~ 1. 

If were we considered only the tangential (^-mode, 
our problem reduces to the case first investigated by 
Sutherland.- - Furthermore, the ip-mode is fully studied 
in the context of the chiral helimagnet i 37 ' 39 However, to 
realize the longitudinal magnetic current, as we will see, 
it is essential to include into consideration the 0-modc. 
Even of zero-field, k = 0, the #-mode acquires the en- 
ergy gap (JS 2 ) {m/n)q = DS 2 ^ The 6>-gap directly 
originates from the Dzyaloshinskii interaction that plays 
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FIG. 3: (a) Spatial variatoin of the gap function A (x) for the 
i9-mode. (b) The minimum gap A m i n is shown as a function 
of H /H c . 



a role of easy plane anisotropy. On the other hand, the 
cp-mode is the massless Goldstone mode corresponding to 
rigid rotation of the whole helix around the helical axis. 40 
Even after switching the perpendicular field, the #-mode 
(t^-mode) remains to be massive (massless). 
The mode expansion is 

v(x) = ^i] a v a (x), u(x) = ^£, a u a (x), (11) 

a a 

where the orthonormal basis v a (x) and u a (x) is deter- 
mined through the eigenvalue equations, 

C v v a (x) = p a v a (x), C u u a (x) = \ a u a (x), (12) 

with a mode index a. The vibrational part is now given 
by 

V = / dx(H u + %,) = V {paVl + ^C) • (13) 
In explicit form the eigensystem (|12[) present the 
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Schrodinger-type equations, 

d 2 v a (x) 



dx 2 



= [2k 2 su 2 (x, k) - (k 2 + 2p a )]v a (x), (14) 



d 2 u a (x) 



- [2re 2 sn 2 (x,k) - (« 2 -4g + 4 + 2X a )]u a (x). 

(15) 



d:v 



In Eq. p5|) we consider the case of weak field correspond- 
ing to small k that admit dn (x, k) ~ 1. In appendix 
A, we present the general scheme to treat the periodic 
potential having the spatial period 2K and show that 
this approximation does not affect qualitative result pre- 
sented below. Now, both equations (TT5|) and (jT5J) reduce 
to the Jacobi form of the Lame equation^ and their so- 
lutions have been discussed by us previously^ 2 - (see also 
Appendix B). The analysis shows that both the <f and d 
mode consist of two bands? 3 -^ i.e., 



Acoustic if mode 



,(-) 



Optical <f mode : 



(16) 



,(+) 



Pa 



(+) 



1 



\/2|sn (a,«')|' 



Acoustic d mode : 



,(-) 



Optical •& mode 

, ,(+) 



2q - 2 



(a,K r ), 



= VAi +) = 



2q - 2 



2sn 2 (a, k') ' 



where the real parameter a runs over K' < a < K'. 
Here, K' means the complete elliptic integral of the first 
kind with the complementary modulus k' = yl — k 2 . 

By imposing the periodic boundary condition, the 
quasi-momentum (Floquet index) is introduced for the 
acoustic 



)(-) 



2KK' 



Z(a,n') , 



(18) 



(0 < \Qi } \ < n/2K) and the optic 



(+) - 



71-0 „ , /x , / is cn (a, k') , „ N 
- + Z (a, k ) + dn (a, k ) 77777^-777 (19) 



sn (a, ft') 



(n/2K < IQ^l) branches, respectively, where Z de- 
notes the Jacobi's Zeta-function. 41 The representation 
was given by Izyumov and Laptev,— and differs from 
a conventional representationi 38 i 41 

A dispersion relation is given by lo as a function of 
Floquet index Q. We show the excitation spectra lu v and 
u)$ in Fig.|U Because 4/ir < q < 2, the energy gap of the 
d mode, A^(a = 0) = ^2q — 2, has a range a/8/7t — 2 < 
A < y2. The gap has a maximum value = \[2 at zero 
field (k — 0) and monotonously decreases as the field 



increases up to the critical field (k — 1). The normalized 
wave function at the bottom of the acoustic band is 



A a=0 (x) 



1 K(k) 1 K(k) 

EML dn {x ' K) =2\ EML d ^ {x) 



(20) 

In the next section we demonstrate that A a= o(a;) exactly 
corresponds to the zero translational mode. 
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FIG. 4: The energy dispersions of the eigenmodes for (a) 
the tangential ^-fluctuation (u} v ) and (b) the longitudinal 8- 
fluctuations 



IV. GALILEAN BOOST OF THE KINK 
CRYSTAL 

In the previous section, we determined the phonon 
modes. Next we consider the translational mode. The 
translational symmetry holds after the kink formation 
and gives rise to the Goldstone mode, i.e. zero mode 
dxfo (x) — 2dn [x, k). Although the Gaussian fluctua- 
tions around the kink crystal state are assumed to be 
small, this is not true for the zero mode which describes 
fluctuations without damping. Then, the center of mass 
coordinate is elevated to the status of the dynamical vari- 
able X(t) and the phonon modes are orthogonal to the 
zero mode. To describe this situation, we follow the col- 
lective coordinate method. 42 ' 43 

At first, we construct the Lagrangian for the kink crys- 
tal system. We make use of the coherent states of spins, 



\m) =exp[*0 i A-S]|S, 5), 



(21) 
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where 



A = 



n x n. t 
I no x n* 



(22) 



with rij = (cos 6i, sinOi cos <pi, sinOi sin ipi) and n = 
(1, 0, 0). Si(i = x, y, z) are the generators of SU(2) in the 
spin-S" representation and satisfy [S a , Sp] = ie a p 1 S 1 . The 
highest weight state \S,S) satisfies Sf \S,S) = S\S,S) 
and S 2 \S, S) = S(S + 1) \S, S). The states |n,) form an 
overcomplete set and give (rij|S|nj) = Srij. Using this 
representation, the Berry phase contribution to the real- 
time Lagrangian per unit area is written as 

£Berry = hS ^(cos 0j - l)d t ipi 
i 

= hS— dx(cos6 -l)dt<p, (23) 
m Jo 

where we took the continuum limit in the second line. 
Now, we construct the Lagrangian, 

£ = c Q [ dx[cos6 - l)d t <p- ciV, (24) 
Jo 

with the coefficients 

c = hS— , ci = JS 2 — , (25) 
m k 

and expand tp and 9 in the form, 

ip(x,t) = ip (x-X(t)) + Y^ioVa{t)v a (x - X(t)) , 
9 (x, t) = vr/2 + £~ £<* (t)u a (x - X(t)) . 

(26) 

In the expansion of the #-mode, it is not necessary to 
exclude a = 0, since the #-mode does not contain zero 
mode. This description amounts to using the curvilin- 
ear basis, {X, r) a , , in functional space and taking the 
generalized coordinates q\ = X, qi a = rj a , q^a = £«• 
Since the zero mode d s ifo (x) is orthogonal to the phonon 
modes, we have 



i dx 
for a ^ 0. Noting that 



dx - v a (x) = 0, 



(27) 



<p = -qi 



oo 

^ q2aV a , 



and 



1 — cos ( 



q3a^a j 



and plugging these expressions into the Lagrangian 
we obtain 



C 



-co 



a 



gi y j fC a q3a (28) 



^2M a ,pq 3a q2p ) - CiV, (29) 

a,p 



where higher order terms O (g 3 ) are dropped. The over- 
lap coefficients are given by 



<Ja Jo dxV a (3?) , 
tCa = Jq dx d(P Q^ u a (x) , 

M a p = J L dx u a (x) vp (x) . 



(30) 



The Lagrangian (128|) is singular because it does not 
contain any term of the form qtqj, and the rank of 
the Hessian matrix UPC/dqidqA becomes zero. This 
means that there is no primary expressible velocities. 
Therefore we need to construct the Hamiltonian by us- 
ing the Dirac's algorithm for the constrained Hamilto- 
nian systems ! 44 ' 45 The details of the treatment have been 
given in our previous treatment (see, also Appendix C). 
The final result is 



0. 



2ci A Q 



(31) 



which means that only finite amplitude of the -d-mode, 



u(x) 



^aU a (x) , 



(32) 



appears when the collective velocity X is finite. This is 
exactly the manifestation of the ODLRO. In other words, 
the u(i)-field is interpreted as the demagnetization field 
that drives the inertial motion of the kink. Using Eq. (|3"Tj) . 
we reach the final form of the physical Hamiltonian, 

^ = ^E A ^^E^ l2 4 Mi ' (33) 

a a 

where the inertial mass of the kink crystal is introduced 

i. /™y 

Z\ v k ) 



2ci 



2 v^ 



(34) 



The physical Hamiltonian (|33[) describes the inertial mo- 
tion of the kink crystal. 

The linear momentum per unit area carried by the kink 
crystal may be presented in the form^ 2 - P = 2irhSQ + 
MX, where the topological charge 



Q = ^MZ)-M0)] 

Z7T 



(35) 



is introduced. Apparently, the transverse magnetic field 
increases a period of the kink crystal lattice and dimin- 
ishes the topological charge Q and therefore it affects 
only the background linear momentum (see discussion in 
Sec. VI). The physical momentum related with a mass 
transport due to the excitations around the kink crystal 
state is generated by the steady movement. 

The "superfluid magnon current" transferred by the 
^-fluctuations is determined through the definition of the 
accumulated magnon density^ p s in the total magnon 
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density N — gp^S (1 — cos 6) = po + Ps, where the su- 
pcrfluid part p s = — gp^S cos 6 is conjugated with the 
magnon time-even current carried by the ^-fluctuations 

fix) = gp B S^- X 2 £ ^ (*) (36) 

ZCl AfV 

a 

via a continuity equation. 22 The important point is that 
the only massive #-mode can carry the longitudinal 
magnon current as a manifestation of ordering in non- 
equilibrium state, i.e., dynamical off-diagonal long range 
order 4^ 

The net magnetization (magnetic dipole moment) in- 
duced by the movement is 

m(x) ~ —gpBSu(x) 

= -gn B S-^-— - u a(%)- ( 37 ) 

a 

The minus sign means that the net magnetization pro- 
duces a demagnetization field. 



V. QUANTITATIVE ESTIMATES 

To compute the mass M, the spin current j x and the 
magnetic dipole moment m, we consider an array of par- 
allel chains described by the model |T]), where a number 
of chains per unit area is n arca = l/ a o- I n the case of 
the molecular-based chiral magnets, the crystal packing 
is usually loose (do — 10 -9 [m]) and the exchange interac- 
tion is rather weak (J ~ 10[K] ~ 1(T 22 [J]). On the other 
hand, in the case of the inorganic chiral magnets, the 
crystal packing is close (do — 10~ 10 [m]) and the exchange 
interaction is rather strong (J ~ 100 [K] ~ 10 _21 [J]). We 
take these values as just typical parameter choices. The 
strength of the Dzyaloshinskii interaction is ambiguous 
and we simply take qo — D/J = 10~ 2 . 



A. Mass of the kink crystal 

The mass M of the kink crystal is given by Eq. ([34|) . 
Evaluation of the overlap integral JC a is performed in 
appendix D and yields K a = 5 a .oK,o, where 



that after simplification yields 



/Co — 2 1 



' E(n)m L 
K (k) k ao 



Therefore we have 



M = 



if. (™) 2 



(38) 



(39) 



The factor 1/a 2 appears here after the MKS units [m] 
for distances are recovered in Eq. ([33|) . The mass per 
unit area is given by 



4 (™\ 2 J_ 
2ci V k J Aq a\ ' 



Afaroa = "area X M = ( - ) — , (40) 



-W*area — 



2E(k) h 2 L H 2 L 



XqK («) Jag Jag 



The last relationship is reliable in the case of small fields, 
i.e. A ~ 2, and K (k) = E («) ~ tt/2. 

Noting that the period of kink measured in lattice 
units is given by Eq.Q, which turns into £ k ink = 
8K(k)E(k) J 'irqo ~ 2irJ/D for small fields, the mass per 
one kink acquires the form 



A/kink = A/ a 



/ki, 



L 



j h 2 

D Jal' 



(41) 



As a typical example of the molecular-based chiral mag- 
nets, we have 

M kink ~ 10- 9 [g/cm 2 ]. 



For the chain length L/ao = 10 5 , we have the total mass 
M aroa ~ -j^M kink ~10 _4 [g/cm 2 ]. As a typical example 
of the inorganic chiral magnets, we have 



M\ 



kink 



10" 6 [g/cm 2 



For the chain length L/ao — 10 6 , we have the total mass 
A^area— 10~ 2 [g/cm 2 ] . This heavy mass should be com- 
pared with the mass of conventional Bloch wall mass in 
ferromagnets. To make clear the difference, in appendix 
E, we gave a brief summary of this issue. In the present 
case, appearance of the heavy mass is easily understood, 
since the kink crystal consists of a macroscopic array of 
large numbers of local kinks. 



B. Spin current 

As it follows from Eq. (|36f the physical dimension of the 
spin current density is Wb • m 2 /s. Using the results of 
the appendix D the spin current density given by Eq. (|36|) 
transforms into 



c m 1 • 2 JC 

J {x) = 9VbS- X z —uo{x) 

IC\ k ao Aq 



(42) 



The factor l/ao occurs after the MKS units for distances 
are recovered in the continuity equation d x — > aod x and 
in the velocity X — > X/a Q . After simplifications with aid 
of Eqs.©, l[2%]). pi)) , and (J3HJ) we immediately have 



^4EW1 . 2d _ 
Jao nq Q A 

For the case of weak fields corresponding to small k this 
yields 



j*(x)~^X 2 dn (x,k). 
Ja qo 



(44) 



We present a schematic view of an instant distribution 
of spins in the current-carrying state in Fig. [5J In Fig. [51 
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we present a snapshot of the position dependence of the 
current density j x {x) in the weak field limit, given by 
Eq. (|44| . In Fig. El we depicted the cases of the magnetic 
field strengths H/H c = 0.1, 0.5, and H/H c ~ 1. Al- 
though the formula (|44| is valid only for the case of weak 
field limit, but qualitative features are well demonstrated 
by just extrapolating the validity up to H / H c ~ 1. As 
the field strength approaches the critical value, the cur- 
rent density is more and more localized. 

For both the typical molecular-based and inorganic chi- 
ral magnets, we have 

f(x) ~ 0.1 hb X 2 ~ 10~ 24 X 2 [Wb • s]. 

Taking the velocity of order X ~ 10 2 [m/s] we obtain 
finally 

f(x) ~ lO" 20 [Wb-m 2 /s], 
therefore the current through the unit area 

Jarea(S) = f X ™arca ~ 1 [Wb/s] . 



FIG. 5: A schematic view of an instant distribution of spins 
in the current-carrying state. This picture corresponds to the 
case of intermediate field strength. 



rm/ti 



H/Hc = 0.1 




H/H c = 0.5 



FIG. 6: A snapshot of the position dependence of the current 
density j x (x). j x (x) is scaled by its maximum j'^ax = J x {0)- 
We depicted the cases of the magnetic field strengths H/H c = 
0.1, 0.5, and H/H c ~ 1. 



i.e. the relationship j x = —mX holds. By the same 
manner as it was made for the spin current we obtain in 
the case of the small fields 



m{x) ~ — ®^ B — X dn (x, k) . 
Ja q 



(45) 



Therefore, for both the molecular-based and inorganic 
chiral magnets, we have 



m(x) r~j 0.1 [IB X r*j 10 /iB, 



(46) 



i.e. m(x) is of order 10 22 [Wb • m]. The total magnetic 
moment of the chain is 



^ B X [ dn (x,k) dx 
Ja qo Jo 



Ja q 



ttQX. 



(47) 



We here used the relations J dn (x, k) dx — am (x, k) 
and [ipo(x) + t]/2 = sin^ 1 [sn (x,k)] — am (x,k) that 
leads to 



dn (x, k) dx = ^[<po(L) - M0)} = *Q, (48) 



where Q is a topological charge introduced in Eq. (T35l) . 
Noting, 



Q — i/Zkink 



TTq L 



8K (k) E («) a ' 
we have the chain magnetization 

gfi B h ( L 



(49) 



^chair 



- x. 



2Ja,Q \ao 
The total moment per unit volume 



m vo\ — m chain x ^arca x L> — 777 I A ■ 

2Ja \a J 

As a typical example of the molecular-based chiral mag- 
nets, we have 

m vo i - 10~ n X ~ 10~ 9 [Wb • m]. 

As a typical example of the inorganic chiral magnets, we 
have 

m vol ~ 10" 8 A - 10~ 6 [Wb- m]. 



VI. DISCUSSIONS OF RELATED TOPICS 



C. Magnetic dipole moment 

The magnetic dipole moment [Eq. (|5T]) ] induced by the 
motion is given by 

/-x Co m 1 • /C , . 

m( x ) = -g^Bd- X —u (x) , 

zci k ao Ao 



A. Background spin current problem: SU(2) gauge 
invariant formulation 

Heurich, Konig and MacDonald 47 proposed that the 
external magnetic fields generate dissipationless spin cur- 
rents in the ground state of systems with spiral magnetic 
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order. Here, we comment on the relevance of the present 
work to this issue. In our model, the background spin 
current is given by 

jb g = d<p (x)/dx- q oc dn(x) - 2E(n)/iv, (50) 

i.e. there arises the misfit of the kink crystal to the 
helimagnetic modulation and consequently the current 
comes up. Below we prove that this current exists on 
a link between two sites but it causes no accumulation 
of magnon density ("magnetic charge") at the site due 
to continuity equation, i.e. the current is not related 
to the magnon transport. This supports reasonings of 
arguments by Schiitz, Kopietz, and M. Kollar— that ap- 
pearance of finite spin currents is direct manifestation of 
quantum correlations in the system, and in the classical 
ground state the spin currents vanish. 

The background spin current problem is best described 
by the SU(2) gauge invariant formulation developed by 
Chandra, Coleman and Larkin. 49 By imposing the local 
SU(2) gauge invariance of the theory, we obtain the fic- 
titious SU(2) gauge fields a and h that give the spin cur- 
rent j( s ) = dCg/da, and the spin density S =dC g /dh, 
respectively, where C g is the gauge-invariant Lagrangian. 

Following Chandra, Coleman and Larkin, we use the 
SU(2) Schwinger boson representation, 



The SU(2) rotation gi gives the rotation of the spin vec- 
tor, 

= exp (-0, • i) S, ~ S, - 4 x S h (56) 

where {I^) v \ — (/•*, v, \ — x, y, z) is the adjoint rep- 
resentation of the Lie algebra of the SO(3) group charac- 
terized by [I^Iu] = £^xh- 

Rewriting the Lagrangian in the gauge invariant form, 
there appears a term 

b'L {9id T gr l ) b' ia = ie 0i * d T ® t ■ S',, (57) 

that leads to introducing the gauge field transformed 
as 



h' 



8,1 



(58) 



2S, (a =1,2) (51) 



where r = it. Introducing the gauge covariant time 
derivative, T> t = dt — hx,we have h 4 = + T> t ®i.The 
fictitious magnetic field Vt©j is induced by the time- 
dependent rotation of the spin reference frame. 

The exchange terms are regrouped in a gauge- invariant 
form, 



where a — (o~ x , o~ y , o~ z ) are the Pauli matrices. In the 
path-integral prescription, the partition function is rep- 
resented as 



2 = J T>b\ a Vb ta V\ exp {- 



C(t)gIt , (52) 



where the Lagrangian is given by 

i 

+ H[S{b\b)] , 



2S* 



(53) 



where TL is the Hamiltonian |T]) written in terms of 
the Schwinger bosons, and r represents the imaginary 
time. The Lagrange multiplier provides the local con- 
straint. The local SU(2) gauge transformation acting on 
the SU(2) doublet, bf = (bf v b£) is given by 



bt~9i\ 



K = gA, 



where 



i(t) = exp 



(54) 



(55) 



J Si • Sj + D • Y s > x s j 



<i,3> 



<h3> 



—J 2_, ex P 

<i,3> 



a • dr I I. 



Si, (59) 



where J =V 'J 2 + D 2 , Xi represents the position of the 
i-th site, and the spin vector potential is introduced as 
a x — [DJ J) e x , a y = a z = 0, corresponding to the model 
([T]). The form of Eg. (fB"9")) indicates that the tangential 
phase angle tpi can be gauged away by the local rotation 
of the angle (D/J)R X i around the x axis. The gauge 
field a is transformed as 



a i ^a! i =e @iI ( a -d Xi & i ), 



(60) 



or a^ = a — V Xi 0; via the gauge covariant space deriva- 
tive V Xi = d Xi + ax . In addition to the physical gauge 
field, (D/J)e x , there appears the fictitious gauge field, 
V Xi 0£, induced by the spatial rotation of the spin refer- 
ence frame. 

The variation of the partition function under a local 
gauge transformation must be zero 
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5Z = J VblVb la VX t exp (- J C g (r)dr^j - + ■ Sh'^j = 0, (61) 



where <5a- Q = — V Xi <5©i, <5h- = Vt<50j. Consequently, 
one obtains the conservation law 

V,, (dCg/da^) - V t {dCg/dh'i) = 0. (62) 

By definition (<9£ 9 /da'-)| a , =a = J ■ Ms the spin current, 

where the gauge field is fixed by the Dzyaloshinskii vec- 
tor. On the other hand dCg/dh^ — —hi, and we finally 
obtain the continuity equation 

V,, Jf } + V t S, = 0, (63) 

where = j\zl i+1 + J^^. In the explicit form, the 
spin current from the site i to i + 1 is given by 

jSi+l = Jx - ( S * X S *+l) + X i t( D X S '+l) X 
= S 2 J sin (ip i+1 - tfi- ip a ) e x . 

For the long-period incommensurate structure {D/J < 
1) this yields in the continuum limit 

The spin current from the site i — 1 to the site i 

jfL, = M (Si x Si_i) - [(D x Si) x Si_i] 

gives —J i _> i+1 in the continuum limit and compensates 
Thus, the spin current through the i-th site causes 
no accumulation of magnon density at the site, i.e. 
the current is not transport one. The accumulation of 
magnon density means that the local quantization axis 
is wobbling. This wobbling motion, however, contradicts 
the spontaneous symmetry breaking in the ground state. 

B. Spin supercurrent in 3 He 

The moving kink crystal belongs to a class of dynam- 
ical systems out of equilibrium. 15 In contrast to a class 
of equilibrium macroscopic ordered state with a broken 
symmetry (ordered magnets, liquid crystals, superfluids 
and superconductors) an emerging steady state is sup- 
ported by pumping of energy. The coherent spin pre- 
cession discovered in superfiuid 3 He known as homoge- 
neously precessing domains (HPD) is a striking example 
of the quantum stated 

The precession of magnetization (spin) occurs after 
the magnetization is deflected by a finite angle by the 
rf field from its equilibrium value. The Larmor preces- 
sion spontaneously acquires a coherent phase throughout 
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the whole sample. This is equivalent to the appearance 
of a coherent superfiuid Bose condensate, i.e. HPD is the 
Bose-condensate of magnons. According to the analogy 
the deviation of the spin projection from its equilibrium 
value in the precession plays the role of the number den- 
sity of magnons. In terms of magnon condensation the 
precession can be viewed as the off-diagonal long-range 
order for magnons, where the phase of precession plays 
the role of the phase of the superfiuid order parameter, 
and the precession frequency plays the role of chemical 
potential. 

The remarkable property of the magnon Bose conden- 
sate in 3 He-B is that non-equilibrium precession has a 
fixed density of Bose condensate. The density cannot re- 
lax continuously, a decay of the condensate occurs due 
to decreasing volume of the superfiuid part. This results 
in the formation of two regions of precession: the do- 
main with HPD is separated by a phase boundary, where 
a precession frequency equals to the Larmor frequency, 
from the domain with static equilibrium magnetization 
(non-precessing domain, NPD). In the absence of a con- 
tinuous pumping, i.e. rf field, HPD remains in the fully 
coherent Bose condensate state, while the phase bound- 
ary between HPD and NPD slowly moves up to decrease 
a volume of the Bose condensate. 

We may suggest that in the total analogy with the 
supercurrents in 3 He, i.e. spin currents transferred by 
the coherent spin precession, the pumping of magnons in 
the kink crystal (by ultrasound, for example) will cause 
an appearance of homogeneously moving domains with 
ODLRO separated by a phase boundary from the domain 
with a static soliton lattice. Without an external flux of 
energy, the relaxation will occurs via gradual decrease of 
the volume of the superfiuid phase. 



C. Experimental aspects 

In realizing the bulk magnetic current proposed here, 
a single crystal of chiral magnets serves as spintronics 
device. The mechanism involves no spin-orbit coupling 
and the effect is not hindered by dephasing. Finally, we 
propose possible experimental methods to trigger off the 
spin current considered here. 



1. Spin torque mechanism and spin current amplification 

The spin-polarized electric current can exert torque to 
ferromagnetic moments through direct transfer of spin 
angular momentum.-. This effect, related with Aharonov- 
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Stern effect^ for a classical motion of magnetic moment in 
an inhomogeneous magnetic field, is eligible to excite the 
sliding motion of the kink crystal by injecting the spin- 
polarized current (polarized electron beam) in the direc- 
tion cither perpendicular or oblique to the chiral axis. 
The spin current transported by the soliton lattice may 
amplify the spin current of the injected carriers. 



2. XMCD 

To detect the magnetic dipole moment dynamically in- 
duced by the kink crystal motion, x-ray magnetic circular 
dichroism (XMCD) may be used. Photon angular mo- 
mentum may be aligned either parallel or anti-parallel to 
the direction of the longitudinal net magnetization. 



by 



EM = 



dt 



3mo 9^b ft 



p.r 



4?r Ja Q q ( x 2 + p 2) 



,2\5/2 ' 



(66) 



where p is the radial coordinate. Then, in the metallic 
chiral magnets, strong energy loss may occur due to the 
induced eddy currents. This phenomena is exactly analo- 
gous to a well known fact that a magnet moving through 
inside of the metallic pipe feels strong friction. On the 
other hand, in the insulating chiral magnets, there is no 
eddy current loss and instead the polariton excitations 
are expected to occur. Therefore, the frictional force act- 
ing on the moving crystal can be strongly diminished in 
the insulating magnets^ 



3. Ultrasound attenuation under the magnetic field 

Further possibility to control and detect the spin cur- 
rent is using a coupling between spins and chiral torsion. 
Fedorov et al^- first pointed out that under the exter- 
nal torsion, the magneto-elastic coupling of the form, 
J^Ri Rj 9ij [ v x (uj — Uj)] ■ Si x Sj, appears, where is 
the displacement of the magnetic atom at a lattice point 
Hi. Then, the quantity dy = [V x (u, — Uj)] plays 
a role of an effective Dzyaloshinskii interaction. Ultra- 
sound with the wavelength being adjusted to the period 
of the kink crystal may resonantly modulate and may 
exert the periodic torque on the kink crystal. Conse- 
quently, the kinetic energy is supplied to the kink crystal 
and the ultrasound attenuation may occur^i Then, the 
attenuation rate should change upon changing the ap- 
plied magnetic field strength. 



4- TOF technique 

The most direct way of detecting the traveling magnon 
density may be winding a sample by a pick-up coil and 
performing the time-of-flight (TOF) experiment. Then, 
the coil should detect a periodic signal induced by the 
magnetic current. 



5. Energy loss of the moving kink crystals 

The moving kink crystal produces the time-varying 
vector potential per kink, 

47T r 6 

where V = X and r is the position vector with respect 
to the kink center. Then, the magnitude of the induced 
azimuthal electric field E around the chiral axis is given 



VII. CONCLUDING REMARKS 

In this paper, we gave a detailed account of a mecha- 
nism of possible longitudinal transport spin current in the 
chiral helimagnet under transverse magnetic field. The 
most important notion is that the "spin phase" directly 
comes up in the observable effects through the soliton lat- 
tice formation. In our mechanism, the current is carried 
by the moving magnetic kink crystal, where the linear 
momentum has a form, P = 2irSQ + MX. The topolog- 
ical magnetic charge, SQ, merely enters the equilibrium 
background momentum 2irSQ, while the collective trans- 
lation of the kinks with the velocity X gives the mass M. 
Among the Gaussian fluctuations around the kink crys- 
tal state in the soliton sector the longitudinal (along with 
the helical axis) fluctuations play a crucial role to deter- 
mine the mass of kinks. Appearance of the spin currents 
is a manifestation of ordering in non-equilibrium state, 
i.e., dynamical off-diagonal long range order. 

We also stressed that if we took account of only the ip- 
fluctuations, the spin current (Josephson current) would 
cause no accumulation of magnon density and the current 
is not transport one. The accumulation of magnon den- 
sity means that the local quantization axis is wobbling 
but this contradicts the spontaneous symmetry breaking 
in the ground state. 

This mechanism is quite analogous to the Doring- 
Becker-Kittel mechanism of the domain wall motion, i.e., 
the Galilean boost of the solitonic kink. In our case, the 
coherent motion of the kink crystal is dynamically in- 
duced by spontaneous emergence of the demagnetization 
field. To describe the kink crystal motion and resultant 
emergence of the demagnetization field, we revisited the 
Sutherland's seminal work2& and generalized it to the 
case of vectorial degrees of freedom, i.e., not only the 
tangential tp but and the longitudinal degrees of free- 
dom are considered. To clarify the physical meaning of 
the inertial mass, we used the canonical formulation of 
the kink crystal motion. We showed that in the case 
of molecular-based chiral magnets, the inertial mass per 
kink amounts to Mkmk — 10 _9 [g/cm 2 ] and the total mass 
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M arca ~ 10 4 [g/cm 2 ]. In the case of the inorganic chi- 
ral magnets, Mkjnk — 10~ 6 [g/cm 2 ], and the total mass 
Af arca ~ 10 _2 [g/cm 2 ]. Furthermore, the magnetic dipole 
moment per kink, induced by the kink crystal motion, 
amounts to m ~ 0.1 fiB A ~ 10 \ib- Appearance of the 
heavy mass is a consequence of the fact that the kink 
crystal consists of a macroscopic array of large numbers 
of local kinks. 

We here mention that in our scheme, the energy gap 
of the 0-mode plays a role of "protector" of the rigid 
sliding motion of the kink crystal. To excite the 0-mode, 
we need to supply the energy via the external force. This 
situation is reminiscent of the existence of a threshold like 
Larmor frequency in the superfluid 3 He. To make clear 
the physical nature of the edge velocity in our scheme 
is beyond the scope of the present work. We leave this 
problem for future consideration. 

Detection of these observable quantities may be quite 
a promising challenge for experimentalists. Behind this 
issue, there is an actively argued problems on how to 
make use of the indirect couplings among the magnetic, 
electronic, and elastic degrees of freedom. For exam- 
ple, magnetic-field-dependent ultrasonic attenuation may 
give us a new insights. To materialize the theoretical 
model presented here, symmetry-adapted material syn- 
thesis would be required. So far, a novel category of 
materials suitable for chiral magnets has been success- 
fully fabricated on purpose for application in the field of 
both molecule-based and inorganic magnetic materials. 
The interplay of crystallographic and magnetic chirality 
plays a key role there. The materials of this category 
are not only of keen scientific interest, but they may also 
open a possible new window for new device synthesis and 
fabrication in spintronics. 
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and Q is a Floquet index. It may be shown (see Ref^ 3 -, for 
example) that boundary points of bands are determined 
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cos(2KQ) = ±1, 
that produces boundary points of Brillouin zones 
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The periodicity condition (|A2[) means that the potential 
may be expanded into the Fourier series 



lG n X 



V(X) = J2 V nC 



where the reciprocal lattice points are G n = 2irn/(2K), 
n is integer. 

To find Fourier coefficients of the potential V{x) we 
use Fourier series for dn (x) and sn 2 (ai) functions, 

. , . 7r 7r cos (imx/K) 

dn(x) = I > - - — 

y ' 2K K t— 1 cosh (irnK'/K) ' 

n— 1 v ' ' 



and 



o, „ K — E n 2 n cos (imx/K) 



Kk 2 ^ k 2 K 2 sinh (imK' IK) ' 

n— 1 v ' ' 

Plugging these series into (|A3[) we obtain 



,2 £ 2i_ 
V = 1 + k -6— + — g , 
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APPENDIX A: PERIODIC POTENTIAL AND 
BLOCH THEOREM 



We have Schrodinger equation 
dx 2 

where the periodic potential has a period 2K 

V(x + 2K) = V(x), 
and given explicitly by 



2 + V{x)u = eu, 



V(x) = 6K 2 sn 2 (x) - k 2 - 4 + 4q dnx. 



(Al) 

(A2) 
(A3) 



V n = - 



3tt 2 



2irq 



K 2 sinh (imK'/K) K cosh (irnK'/K) ' 



The zeroth-order component Vq determines a shift and 
may be omitted while the component V n mixes the plane 
waves with wave vectors k and k' = k + Tm/K 



k'\v(x)\k) = Vnh,,u +Gn - E v nh'u™ 



IK- 



Hence, a quasidegenerate perturbation theory built in the 
subspace spanned by two states \k) and \k + G n ) 



El Vn 

V* 



= 0, 



(A4) 



yields bands 



k k+G n 
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The gap between the states 
is 



■Qfel) and | -Q^l + Gr, 



2IK 



6tt 2 



4Trq 



1 



K 2 sinh (nnK'/K) K cosh (irnK'/K) 
and it falls rapidly to zero with increasing of n 
2\V n \ ~ exp (-imK'/K). 

APPENDIX B: LAME EQUATION 



The basic properties of the Lame equation are pre- 
sented here. We start with the Jacobi form which is 
defined by 41 



d 2 Aa{x) 
dx 2 



+ l) K 2 sn 2 (x, k) - K + 2 (1 + A)] A & (x), 

• » * 1 • 

where I = 1 and A being a constant. The spectrum is 
labeled by a complex parameter a and given by 



Ar 



1 



r dn a. 



(B2) 



The solution of the Lame equation is exactly given in the 
conventional form 41 - 



_ H(£-a) z(e) 
9(.t) 



(B3) 



where H, 0, and Z are Jacobi's eta, theta, and zeta func- 
tions, respectively, with the elliptic modulus k. Now, we 
require Eq. (IB3|) to be a propagating Bloch wave, i.e., 
Z{a) to be pure imaginary. Recalling that the zeta func- 
tion Z{a) is singly periodic with the period 2K, we see 
that two segments (K - 2iK',K] and [-2iK',0) for a 
are sufficient to fully describe the solution (|B3j) . 
Because of the quasi-periodicity, 

H(x + 2K-a) = -if (a; - a), Q(x + 2K) = Q(x), 

we have 

A & (x + 2K) = -e 2KZ ^A a {x), 
and it is convenient to introduce the Floquet index 



Q{a) = — + iZ(a,k). 



Then, we have 



A a (x + 2K) = e~ 2KlQ A a (x), 

that is analogous to the Bloch theorem where 2K and Q 
have the meanings of the lattice constant and the quasi- 
momentum, respectively. Furthermore, imposing the pe- 
riodic boundary condition 



Aoc{x) 
(B4) 





' e -2KiQ 







A a (x + L) = A a la; 

= e- iLQ A 6l (x)=A & {x), 



we have the quasi-momentum as usual, 
Q = -jr-n, 

where n is integer. 

Finally, we have the Bloch form, 



A a (x) 



HQe a ) r -iQx i^x 

Q(x) 



Other than the conventional parameterization, it is con- 
venient to work with a real parameter a related with a 
by 



a — ia + K — iK' , 
for the acoustic branch, and 

a — ia — iK' , 



(B5) 
(B6) 



for the optic one. Within the new parametrization the 
eigenfunction for the acoustic mode transforms in the 
following way 



H(x -ia-K + iK') 



where xq — ia + K, and (i = 1,2,3,4) denote the 
Theta functions. Furthermore, we have 



^3 (^N - iK']) = e^e*** ^ (^x ) . 



and 



e*(£[x -iKi) 



R(x -a) = i \ '-er^di ( — 



§2 (£c x ° 



\2K 



X — Xq 



that yields 



A a (x) = i 



^(^[xo-iK'])^(^[x-x [ 



(B7) 



This is an alternative representation for the solution 
(|B3j) . 37 and it is used in the paper. The case of the optic 
branch (xq = ia) is considered by a similar way. 

The transformation of the Floquet index for the acous- 
tic branch is carried out as follows. By noticing that 

Z(ia + K -iK') 

= Z{ia) + Z(K - iK') - sn (ia)dc(ia) 

. . KE' + K'E - KK' , . , , . 

= Z(iot) + i — sn umdcua) 

K 

= i iZia. k') — iir — — — , 

2K v ' ' 2KK' 1 
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where we used the Jacobi's imaginary transformations 
and the Legendre's relation KE' + K'E - KK' = it/2. 
Therefore, we have 



Q(a) = Q(a) 



2KK 1 



Z(a, n'). 



The same transformation for the optic mode (a = ia 
iK 1 ) yields 



Q(a) 



2KK' 



Z (a, k ) + an (a, k ) ; — . (tie,) 

sn (a, k 1 ) 



By the same manner, the corresponding spectrum is 
parametrized as 



A a = — -dn a = A a = 



^5-sn 2 a (acoustic) 
1 

~2 T~ (°P tlc ) 



Now, we briefly review the origin of the band 
structured In the limit k — > 1, the Lame equation re- 
duces to the Schrodinger equation, 



where E 
potential 



d 2 ip(x) 
dx 2 

■■ k 2 (l + 



f E + U a sech 2 (ax) = 0, 
A) - 1(1 + 1), U Q = £( 



l).The 



U(x) = — U sech 2 (ax), 



is modified Poschl- Teller potential and for I = 1, there 
are one bound state and one 'perfectly transmitted (re- 
flectiordess) scattering stated The band structure of 
the Lame equation is understood as follows. In the limit 
of well separated modified Poschl- Teller potential, the I 
bound states give discrete levels and the scattering states 
give broad continuum. When the potentials form a lat- 
tice, the discrete level overlaps and the energy band may 
be formed. Even after the band formation, the gap be- 
tween the bound level and the scattering continuum re- 
tains. Therefore, the resulting band is split into the lower 
acoustic band and the upper optical band. 



APPENDIX C: DIRAC'S CANONICAL 
FORMULATION FOR THE SINGULAR 
LAGRANGIAN THEORY 



The Lagrangian (|28[) itself gives rise to a set of primary 
constraints, 




= Pi - CO Ta £»q3a ~ 0, 
a = Via + C (j n + Tfi M a (3q3f3 
- P3a « 0, 



(C3) 



where the symbol ps means "weakly zero," i.e. 
may have nonvanishing canonical Poisson brackets with 
some canonical variables. Because of a lack of primary 
expressible velocities the Hamiltonian with the imposed 
constraints, 



H* 



?2c 



E^&a+ClV, (C4) 



coincides with H c , i.e. q\, q2 a , and q^ a (primary inex- 
pressible velocities) play the role of Lagrangian multipli- 
ers. Now, the Hamiltonian H* governs the equations of 
motion of the constrained system. The relevant non-zero 
Poisson brackets are computed as, 



{4a. 03a} = CoTfiM a fi, 

= -2cip a q 2a , 

= — 2ciAa93a, 




(C5) 



and {qi,pj} = Sij gives rise to the constraint conditions, 
= {^\H*} = 0, $1 = {4a> ff *} = °. and 
03^ = {42) H*^ = 0, or in the explicit form 



C0^]^a?3Q = 0, (C6) 

a 

Co E Mapqzp - 2c 1 p a q 2a = 0, (C7) 





The canonical momenta conjugate to the coordinates 
X(t), r) a (t), and £ Q (t) are given by 

pi = dC/dqi = c T a fcaq3u, 
p 2a = dC/dq 2a = -c (j a + TpM a pq3fjj , (CI) 

p 3a = d£/dq3a = 0, 

and we obtain a canonical Hamiltonian, 

H c = Piqi + ^P2a92a + ^ C. (C2) 



\5o< 



Co K-a<ii - ^2,M a fiq2a - 2ciA 

aQ3a 



= 0. (C8) 
Eq. (|C6|) gives q 3a = and then Eq. (|C7|) gives q2 a = 0. 

(2) 

Now, there arises the secondary constraints (jr a = q2 a ~ 

•(2) ( (2) 1 

to be constant in time, <p a ' = < <f> a , H* > = q2 a = 0, 

and the consistency condition is fulfilled. Finally Eq. 
(|C8[) relates q^ a = £ Q to q\ = X, and produces Eq.d3TJ). 
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APPENDIX D: COMPUTATION OF K, a 



We compute 



JC a = 2 / dx dn (x, n) u a (x) 
Jo 

= 2N(a) dxdn(x, K ) KZK n ' 

with xq — ia + K, and where N (a) is a normalization 

tactor. Noting that dn (x, k) — ^ * — - has a period 

2_fiT, we perform the Fourier decomposition, 

#4 (^(^-^o)) 



dn (a;, k) 



E 



7/e 



where the coefficients arc evaluated as 



1 



K 



7; = / ofedn (x, k) 

2KJ_ k 



Then, we have 



y2K 



K a = 2N(a)Lj2jlh,%i- 
i 

Within the acoustic branch (0 < \Q\ < ^ ), only Q = 
(a = 0) contributes to IC a . Eventually, the orthogonality 
condition ([2T]) of a denumerable basis enforces that there 
is no contribution of the term with I ^ 0. By using 

A a= o(x) — J le(1) ^ n ( x i K )i therefore we have JC a — 
Sa.oK-o, where 



/C = 2* 



' K(k) 
E(n)LJ 



dn 2 (a;, k) dx = 2-t 



K(k) ' 



where we exploited the relation E (k) 
Jq K ^ dn 2 (x, k) dx. 



7Ta/ JS 2 /Ka,Q with if denoting the anisotropy energy. 
Now, let us consider the Bloch wall formed along the 
x-axis and spins are confined to the j/z-plane that winds 
180°. Ddring proposed that the translation of the domain 
wall is driven by the appearance of the local demagne- 
tization field H x inside the wall that violates the condi- 
tion V-M = 0, i.e. H x = —4ir [M x — M x (oo)], and causes 
the precessional motion of the magnetization within the 
yz-plane. Then, the corresponding Larmor frequency 
amounts to cjl = <p = jH x , where 7 is a gyromagnetic 
ratio. On the other hand, in the steady movement of 
the wall, ip = — (d x ip) V, with V being the velocity, and 
consequently we have 

H x = -7- 1 (d x <p) V. (E2) 

The excess of magnetization energy 



gives the energy stored in the moving wall. Taking the 
form AW — AfD6ring^ 2 /2 the inertial mass of the wall 
first proposed by Doringi 7 - is introduced 



A^Doring 



4?T7 2 



d(p 
dx 



dx. (E3) 



The explicit form (dip/dx) depends on kind of domain 
walls, their orientation around crystallographic axes, for 
example 



Mi 



1 



Doring 



47T7' 



■VK/J 



for 180-degree domain wall parallel to crystallographic 
plane (100). Taking into account that 7 = 1.84 x 10 7 
[(Oe • s) _1 ] this yields in the case of Fe 

Mnermg =s 10- 10 [g/cm 2 ]. 



APPENDIX E: 



INERTIAL MOTION OF BLOCH 
WALL 



We here discuss the relevance of the present for- 
mulation to the Doring-Becker-Kittel mechanisim 17 i 18 i 19 
We consider a conventional Bloch wall in ferromagnets, 
where the magnetization rotates through the plane of the 
wall. The wall size is determined by the exchange energy 
cost and the anisotropy energy that amount to 



K 2 JS 2 



KNa, 



o- 



(El) 



where TV is the number of spins inside the wall. Min- 
imizing this energy leads to the wall size Zeioch = 



From Eq. (fE2J) it stems that 



M x = - H x 
Air 



provided M^oo) = 0. This equation should be compared 
with Eq. ([37|) rewritten in the form 



m(x) 



gush 

2Ja qo 



We see that mix) may be interpreted as the demagne- 
tization field and the physical Hamiltonian (|33p may be 
regarded as the energy cost associated with the demag- 
netization process. 
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